One Dimensional Phonon Coupled Electron Tunneling: A 

Realistic Model 

E.Pazy and B.Laikhtman 

Racah Institute of Physics, Hebrew University, Jerusalem 91904, Israel 

Abstract 



The transition probability for a one dimensional tunneling electron cou- 
pled to acoustical phonons is calculated, with the Feynman path-integral 
method. We considered a realistic electron phonon interaction (deforma- 
tion potential, piezoelectric), making use of slowness of the phonon system 
compared to electron tunneling. We show that the problem of the complex 
non-linear coupling of a tunneling electron to the zero point fluctuations of a 
phonon field is equivalent to that of an electron tunneling through a slow fluc- 
tuating spatially uniform barrier, thus resulting in an increase of the tunneling 
probability due to electron coupling with zero-point phonon oscillations. We 
calculated also the energy change of the tunneling electron due to phonon 
emission. 



I. INTRODUCTION 

Tunneling, being one of the most remarkable non-classical manifestations of quantum 
dynamics, has been largely investigated. One of the major issues in the field is the under- 
standing of the properties of tunneling particles coupled to degrees of freedom representing 
the environment. The case of a particle or macroscopic degree of freedom tunneling out 
of a metastable state in the presence of arbitrary linear dissipation mechanism, has been 



extensively studied in a series of papers by Caldeira and Leggett |T[] and other groups [||. 
Two state systems coupled to a dissipative environment have been studied in depth by a 
large number of groups [[J, in the last two decades. 

The problem that we address in the present paper is the tunneling of an electron with a 
given energy E across a rectangular barrier where it is coupled with phonons. We calculate 
the change of the tunneling probability due to this coupling and the average energy loss due 
to phonon emission of tunneled electrons. In this work we consider only a one-dimensional 
problem that can be realized as a quantum wire, or in narrow constrictions. 

This problem, finding the transmission probability and spectrum of a tunneling electron 
coupled to phonons is of fundamental interest as well as being potentially of great tech- 
nological importance. Examples are tunneling through Josephson junctions [|J, bit errors 
in mesoscopic logic and memory circuits ||, as well as quantum cascade lasers |J, and of 
course, scanning tunneling microscopy (STM) 0. There are even some recent examples 
of the use of the transmission spectrum to define an experimental test between quantum 
mechanics and a class of alternative theories ||. Also recent development of nanostructure 
technology makes it possible to check experimentally the effect of coupling with the envi- 
ronment on the dependence of the tunneling probability on the width and the height of the 
barrier. 

Typically the coupling between the tunneling particle and the environment is assumed to 
be linear with respect to particle coordinate. This assumption may be justified for tunneling 
of heavy atoms or macroscopic degrees of freedom but it is not realistic for electrons. A non- 
linear coupling with phonons was considered by Bruinsma and Bak [|9] . But they considered 
only the case when the coupling between the tunneling particle and each of phonon degree 
of freedom is the same, which is not realistic in a regular lattice. Due to this the result of 
Bruinsma and Bak contains a cut off frequency. In our realistic treatment of electron-phonon 
coupling a characteristic frequency naturally arises in the theory. Since the only length scale 
in the problem is the barrier length and from it a typical time scale emerges. 

A recent treatment of the problem was given by M. Ueda [10|, who used the effective 



action technique and got analytic results, for some cases of the spectral density function 
of the phonon bath. Whereas Bruinsma and Bak |9j calculated the transmission spectrum, 
we show that the main effect, on the electron transition probability, is static i.e a static 
effective lowering of the barrier by phonons. Such an effect, a static lowering of the barrier 
due to coupling of the tunneling electron to zero point fluctuations, was also discussed in 
two papers by Ueda and Ando [11J although they studied a different physical problem, the 



coupling of the tunneling electron to electromagnetic modes (in their paper this was referred 
to as the dynamic limit). 

A deformation potential coupling of tunneling electrons with phonons was also studied 
by Gelfand et al. Jl2|]. However, they considered only an infinitesimally thin barrier. 

We consider interaction of tunneling electron with acoustic phonons via realistic cou- 
pling mechanisms: deformation potential or piezoelectric. Coupling with optical phonons 
is neglected for the following reasons. The amplitude for zero point fluctuations for optical 
phonons, for a given frequency is proportional to one over the square root of the frequency. 
Since the optical frequencies are large compared to those of acoustic phonons we consider 
their zero point fluctuations as negligible. Also the coupling constant with optical phonons 
in III-V materials [TB|] is not very large giving a further justification for neglecting the effect 



of optical phonons. Concerning optical phonon emission and absorption we study the low 
temperature case and assume that electron energy and temperature are much smaller than 
the energy of the optical phonons so that emission and absorption can be neglected. 

The coupling with acoustic phonons leads to a modulation of the barrier. This effect 
exists even at zero temperature when the modulation is controlled by zero-point phonon 
vibrations. In other words, as a result of coupling with phonons the height of the potential 
barrier is not any more a classical parameter determined by fabrication of the structure but a 
quantum mechanical dynamical variable that is characterized by a wave function. Thus the 
complex non-linear problem of a tunneling electron coupled to a phonon field is transformed 
to a that of an electron tunneling through an effective potential barrier where the barrier 
height fluctuates. 



The paper is composed in the following way. In Secjlj we formulate the model. In 
Sec.jTTT] we introduce two main approximations which we use to facilitate the calculation of 
the tunneling probability. Both coupling mechanisms that we consider are nonlinear with 
respect to the electron coordinate. This nonlinearity makes the application of generally used 
instanton method difficult and to solve the problem we developed an expansion on some 
small parameters. First, it appears that frequencies of typical phonons interacting with a 
tunneling electron are so small that the phonon system can be considered slow compared 
to the electron. Second, we use the WKB approximation. The slowness of the phonon 
system allows us in the first approximation to neglect phonon dynamics and to consider their 
potential as time independent that dramatically simplifies the calculation (the calculation 
without this approximation appears in the Appendix. The corresponding calculation of 
the transmission coefficient are carried out in Sec.|TV[ The static approximation, however, 
does not include phonon emission and resulted energy loss. To calculate it a more general 
approach is necessary. In Sec.|V] again making use of the slowness of the phonon system 
we go beyond the static approximation and calculate the average energy loss of tunneled 
electrons. 



II. GENERAL FORMULATION OF THE PROBLEM 

A. The model 

We treat the problem of a one dimensional electron coupled to acoustic phonons while 
tunneling through a rectangular barrier, of length L. The potentials restricting the electron's 
movement do not effect phonons, which move freely through the bulk and thus are treated as 
three dimensional. Being interested in the effect of the zero point fluctuations of the phonon 
field on the electron tunneling probability we consider the system at zero temperature. 
Practically it means that temperature is low enough and this condition as well as the finite 
temperature case are discussed in the Appendix. Coupling the electron to the phonon field 
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transforms the problem from a low dimensional quantum problem to a field theory problem. 
A convenient way to approach such a field theory problem is via the formalism of path 



integrals [T^] 



The first stage in the calculation of the transition probability for the tunneling process is 

the calculation of the retarded Green function. We express the energy dependent retarded 

Green function as the Fourier transform of the time dependent propagator, 

roc (%E-T\ 
K(x b ,Uq f ,E | x a ,w qi ) = J K(x 6 ,w q/ ,T | x a ,w qi )exp ( — - — JdT , (2.1) 

that can be written as the combined path integral of the electron and phonons: 

K (x b ,Uq f ,T | x a ,u n i) = J Vx J £>?i q exp [S a (x) + S ph (uq) + <S int (u q ,x)] j , (2.2) 

where x is the electron coordinate, w q are the Fourier coefficients of u(x,t), which is the x 
component of the displacement vector and x a , -u q j ; x b) w q / are the initial and final coordinates 
of the electron and phonons respectively. The limits over the energy Fourier transform result 
from the fact we are calculating the retarded Green function. S e i is the action of the electron 
in the absence of phonons, given by: 

Sd = jf M (f ^ ~ V ) ' (2 ' 3) 
where, m, is the effective electron mass, and V is the constant potential height of the 
rectangular barrier. S p h is the action of the uncoupled phonon field: 

S ph = £dtY, (f I K I 2 - \ ■ " q 2 \ «, I 2 ) , (2.4) 

where p is the crystal density and w_ q = w q *. Since we are treating acoustic phonons 
uj q = q • s, s being the phonon velocity. Si nt is the action associated with the interaction 
between the electron and the phonon environment: 

S^t = - C dtV int . (2.5) 
Jo 

We will treat two electron phonon coupling mechanisms, piezoelectric and deformation po- 
tential, and consider crystals of cubic symmetry. We consider this symmetry as the most 



important because most of III-V compounds belong to this class [|13|] . So for the piezoelectric 
coupling 

Vfnt Z ° = —/tt= H Eq„u q „M q± exp(iq x x), (2.6) 

V ''vol qi/ 

where 

_ Attc 

"qi/ = PijkQiQj^k (2-7) 

is the piezoelectric modules, e is the electron charge, e is the dielectric constant, e v k is 
the polarization vector of the i/-th phonon branch, and V vo i is the normalization volume. 
M q± is the matrix element of the phonon exponent, exp(zq_i_r_i_), between the wave functions 
describing the electron quantization in the cross section of the wire, and r_|_ are the phonon 
wave vector and electron radius vector in the cross section plane. In cubic crystals the only 
nonzero components of the piezoelectric module are f3 xyz = (3 yxz and all other permutations 
fl"5|j . They are equal and we will designate them as (3. 

The deformation potential couples electrons only with the longitudinal phonon mode and 

Vint = —ffr= Yj «A|g|M qi M q± exp(iq x x) . (2.8) 

V Vvol q 

Here A is the deformation potential constant. 

We neglect the screening of V^ nt . Underneath the potential barrier there are no electrons 
and the screening by remote electrons is small. 

B. Transition amplitude and transition probability 

We can now define the transition amplitude through the use of the retarded Green 



function ( |2.2j ). K(xb,u^f, E | x a ,Uqi) is the amplitude to go from initial phonon coordinate 
Uqi to final coordinate w q j and from electron coordinate x a to x^ , for a given energy E, of the 
joint system. (In all intermediate calculation we suppress the index of the phonon branch to 
simplify the notations.) The transition amplitude for the joint system to go from a phonon 
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state designated by ipo (w qi -) to a state ij) n (w q /) is expressed by K(x b , w q /, E \ x a , w q i) in the 
following manner : 

-4nb,0a = J dUqidUqflf)* (« q/ ) iT(x 6 , W q/ , £ | £ a , M qi )^0 (Wqi) • (2.9) 

We treat the problem with the temperature equal to zero, therefore phonons are assumed 
to be initially in the ground state: ipo (w q i). 

The transition probability is the absolute value squared of the transition matrix element. 
Since we are not interested in the final phonon configuration, final phonon states are summed 
over. From the completeness relation we get the following delta function: 5 (u q f — £t q /). We 
end up with the following expression for the electron transmission probability: 

Pb,a = X] I A n b,0a | 2 — 
n 

= J dUqi du^i dUqf % (Uqf) K(x b , U q /, E I X a U OLi )k*(x b ,U CLf ,E I XaU^lpQ (ttqj) . (2.10) 

This equation defines the transmission coefficient. 

III. MAIN APPROXIMATIONS 

The usual approach to the calculation of the functional integral in Eq. ( |2.2p that has been 
started by Caldeira and Leggett [|lj is the integration out the phonon degrees of freedom. 
This can be done exactly because the action is quadratic in u q . In the realistic electron- 
phonon coupling, this results in a complicated electron-electron effective potential, since 
the coupling mechanism is non-linear in the electron coordinate, thus a simpler approach is 
needed. 

We use the following approach. We consider the situation (which is typical experimen- 
tally) when (i)the barrier is so high that the tunneling can be considered semi-classically and 
(ii)the interaction energy between electrons and phonons is small compared to the height 
of the barrier. The first point allows us to make use of the semi-classical approximation to 
integrate with respect to all electron paths. The main contribution to this integral comes 
from only one optimal trajectory that satisfies the equation 
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dV int 

mx = -— , (3.1) 
(for a rectangular barrier dV/dx = 0). This equation has to be solved for given u(x, t). The 



second point allows us to consider Vi n t in Eq.( |3.1|) as a small perturbation and thus it can 



be neglected. As usual in tunneling problems, p6yi7| we change the time t to —it (and T 
to iT) so that the first integral of Eq. ( |3.1| ) has the form 

^ = V-B. (3.2) 

Here E' is the integration constant that can be considered as the electron energy. The 
value of E' is determined from the boundary conditions x(0) = x a and x(T) = xj,. For a 
rectangular barrier when x a and Xj, are fixed E' is a function of T only. Without phonon 
emission E' eventually appears equal to the energy of the incident electron E that is also the 
total energy. If during tunneling phonons are emitted then the electron energy E' appears 
smaller than the total energy. 

Thus, in the tunneling problem an important parameter appears the electron velocity 



v = J2(V — E') /m. If V — E' is around 10 meV or larger than in GaAs where m = 0.067mo 
(mo is the free electron mass) v is around or larger than 2 x 10 7 cm/s. This velocity is larger 
than the sound velocity s by about two orders of magnitude. That means that frequencies 
of phonons with the wave length around the length of the barrier L are much smaller than 
the inverse time necessary for an electron to traverse the barrier. Actually the wave length 
of a typical phonon interacting with the tunneling electron is limited not by the length 
of the barrier but by the width of the constriction a that is smaller that L. However, 
practically the ratio a/L for a constriction where tunneling is still measurable is not very 
small and the assumption that the typical phonon frequencies are much smaller than the 
inverse traverse time is justified. This assumption means that during the traverse time 
(L/v) the phonon potential nearly does not change. In the extreme case we can consider 
it constant. We call this case the static approximation. A similar approximation was used 
by Flynn and Stoneham and later by Kagan and Klinger [F9] treating the problem of 



quantum diffusion of atomic particles. 
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It should be mentioned that the physical significance of the parameter uto, to being the 



phonon frequency, to the tunneling time as defined by Biittiker and Landauer |20], t = — 
(in our case of a rectangular barrier) was also noted in the two papers of Ueda and Ando 
||1 1|| , but in their problem the typical frequency uj was a tunable parameter defined by the 
properties of the electric circuit. 

In the static approximation the problem is dramatically simplified because we need to 
consider electron motion in a static potential and the stationary Schrodinger equation is 
enough for this. The integration with respect to phonon coordinates is reduced to the 
integration only with respect to u^. The derivation of the corresponding expression for 
the transmission coefficient in the static approximation from Eq. (|2.10| ) with the help of the 
expansion in s/v is given in the Appendix. 

In the static approximation we can add also V^f" = V^ nt (M qi ) to the rhs of Eq.(|372|). That 
immediately shows that the small parameter characterizing the interaction with phonons is 
Vint/{V — E). In the calculation of the transmission coefficient we take into account only 
terms of the first order in this parameter. So we ignore corrections to the trajectory x = v t 
coming from the interaction with phonons. It is obvious in the calculation of Vi n t and this is 
also true for S e i. Indeed, the trajectory x = v t is found from the minimization of S e i and any 
correction to this functional contains a correction to the trajectory squared. Such an ap- 
proximation means, in particular, that we neglect all polaron effects. The physical meaning 
of the phonon effect in this approximation is that different phonon configurations change 
the barrier and the main contribution to the transmission comes from the configurations 
corresponding to the barrier being a bit lower in average, so that the tunneling probability 
is higher. It is worthwhile to note that the average (over configurations) height of the barrier 
does not change because (V int ) = but nevertheless the average of the tunneling exponent 
is modified due to V int (similar to, e.g., (exp{V int )) > 1). 

The tunneling across a static barrier is an elastic process and the energy of the incident 
electron and the tunneled one is the same. Phonon emission and the corresponding change 
of the electron energy come about only in the first approximation in s/v when the phonon 



dynamics is taken into account. The calculation of the dynamic correction to the trans- 
mission coefficient that requires a more complete treatment of the modes presented in the 
previous section is given in Sec.|V|. 



IV. STATIC APPROXIMATION 

In the static approximation the problem of tunneling can be formulated in a very simple 
way, without making use of Eq. (|2.2j ). First we can find the transmission probability for a 
given static phonon field. The regular WKB approximation gives for it 

D(E; u qi ) = ^ A{ y + - E f ^ {~l Jo V 2 ™^ + ^ - E) dx} . (4.1) 

The calculation of the electron transmission coefficient in the case when phonon field is 
initially at the ground state is reduced now to the integration of D(E; w qi ) multiplied by the 
ground state phonon wave function squared with respect to « q j, 

P(E) = C 2 J D(E;u^)exp ^-i^ pgs | Mqi \ 2 ^j du^ (4.2) 

(C is a normalizing constant). The expansion in Vi n t gives 

P(E) = P (E)C 2 J expj- ^S= £ V&*dx - f E ? I «* I 2 } *** » ( 4 - 3 ) 

where 

is the transmission coefficient without interaction with phonons. Eq. (|4.3|) can be obtained 
directly from Eq. (|2.10|) (see Appendix). The result can be written in the form, 

P(E) = P (E) exp || $ stat (E)} . (4.5) 

For the deformation potential 

$stat{E) = - f L 9 f |M q jV dV , (4.6) 



8n 2 psiV 
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where si is the longitudinal phonon velocity. 
For the piezoelectric coupling 

^tat(E) -^- 2 J q dq x ^-J |^ q ,| |M q J _ . (4.7) 

Because of the anisotropy the result depends on the tunneling direction with respect to 
crystalographic axes. For the tunneling in [100] the contribution of the longitudinal phonons 
in $ s t a t is small compared to that of transverse ones in a/L and 



= ^ / |M^«M dV , (4.8) 
pv t s t J 



where St is the velocity of the transverse phonons. For the tunneling in [110] direction the 
contributions of both longitudinal and transverse phonons are of the same order, 



= I |M qj 



it + Hi ( l 1 \ Hi 



s t q± \s t sj ql 



rfV (4-9) 



As can be seen the typical phonon wave numbers with which the electron interacts 
are fixed through the length scales of the problem, the barrier length L and the width of 
the constriction a. It is important to note that, the static correction, in the exponent is 
positive, therefore static phonons reduce the effective potential barrier height, enhancing 
electron tunneling probability. The electron, due to coupling to the zero point fluctuations 
of the phonon field, tunnels through an effectively lower potential barrier. 



V. DYNAMIC CORRECTIONS AND ENERGY LOSS 

To get the electron's energy loss due to phonon emission during the tunneling process, 
one needs to go beyond the static approximation. As the first step we make use of the WKB 
approximation to simplify the expression for the propagator K(xb,Uqf,T \ x a) u n i) (|2.2|). In 
this approximation the main contribution to the integral with respect to x(t) comes from 
the saddle point trajectory that in the first approximation in Vi nt /(V — E) is determined by 
Eq. (|3.2|) , As a result the propagator is factorized. In Eq.(|3.2|) we passed to the imaginary 
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time so the trajectory is xo(t) = vt where v = y 2(V — E')/m, thus the propagator can be 
expressed as: 



K(x b ,u cl f,-iT | x aj Uqi) = K (x b ,T \ x a )i^(?Xq/, T 



u. 



qv • 



(5.1) 



Here 
K (x b ,T I x a 



in 



1/4 



exp (— exp f-~ ^ dxyj2m(V - (5.2) 



4 (V + Vgf - E) 
is the electron propagator without interaction with phonons and 

K ph (uqf, T | Uqj) = y X?u q exp j-- [S p h{uq) + «S int (w q , x )] 



(5.3) 



is the phonon part of the propagator. 

So as the electron trajectory is determined the phonon part of the propagator is the 
propagator of an ensemble of forced harmonic oscillators. [|14]] The integration with respect 



to M q leads to 



K P h{uqf, T\uqi) = g{T) exp 



(5.4) 



where 



gen = n 



q y 27r/i sinh tu g T 



(5.5) 



and 



Scl(Uqf,T I M r 



5 I ""qi / 



E 



2 sinh (u^T) 



| COSh(^gT) (| Mqi | 2 + | Mq/ | 2 ) - («qi«q/ + Wqi^q/ 



' 1 rL 



dxfq{x) sinh( 



u 



—u 



If 



\ pU q V JO 
2 rL 



dtf*(x) sinh( 



UJ q X, 

V ' 



1 dxfq(x) sinh(— — 

u q {L - x) 



pUOqV JO 
1 rL 



pUJ q V JO 



dxf*(x) sinh( 



v 



, , dx f X dy \Ux)f*(y) + f*(x)fM] sinh(^) smh( u{L x) 

p z UJq Z V z Jo JO J V V 



)} • (5-6) 



Here according to Eqs. (|2.8|) and (|2.6|) 



IM. 



(5.7) 



vol 
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for the piezoelectric interaction and 

/,(*)= iM^^e**- (5.8) 
v v vo i 

for the deformation potential interaction. One should note that in Eg. ( |5.6| ) we made a 
transformation of variables from t in the action, to x = vt. 

Making use of the factorization ( f5.lD the expression for tunneling probability ( p.lOj ) can 
be written as 

P(E) = J dT J dTK (x b , T | x a )K (x b , f | x a )i ph (T; f) , (5.9) 

where 

£ ph (T;f)=g(T)g{f) 



x J du^ J du^f J du^exp ~ S d (w q /,T | u qi ) - ^ <S d (w q /,T, | u qi 



V'O (Mqi)^o ( U q«) > ( 5 - 



and -00 ( M qi) is the phonon ground state wave function. The integrals with respect to T 
and T in Eq. (|5.9| ) are calculated by the saddle point method [|l£|. Due to the symmetry 
of ^ p h T\ T with respect to the transposition of T and T the saddle point values of these 
variables are identical. We are interested only in the exponential part of the transition 
probability and for this only £ p h [T; T] is necessary. The result can be written in the form 

"2 



i P h (T; T) = A exp 



(5.11) 



where <& p h(T) is calculated in the Appendix and the pre-exponential factor A, will not be 
calculated. 

We now break $ p / l (T) into a static part, Q s t a t, and a dynamical correction so that <& p h = 
&stat + &dyn- As one can expect the static part is identical to $ sta t(E') obtained in Sec.fV], 
Eq.( [4.8|) , in a more simple way. The dynamical correction, &dyn{T, E'), is obtained in the 
Appendix in the leading order in s/v <C 1. Using these notations the saddle point equation 
for time integration is given by: 

— — --^ — -T ^ + £-£' = 0, 5.12 

dE' ) dT dT ' v ; 
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where ® e i{E') = L^2m{V - E'). 

According to Sec.|TV], & e i(E') — $ stat (E') is the electron action in the phonon static field 
at the trajectory with the energy E' . The derivative of the action with respect to the energy 
is the traveling time along the trajectory. The first term, d<& e i/dE', has been defined by 
Biittiker and Landauer as a semiclassical traverse time. The second term, d<& s t a t/dE', 
gives a phonon correction to the traverse time. The sum to the two terms equals T and the 



expression in the parentheses is identically zero. Then Eq.( 5.12 ) can be written in the form 

E-E' = ( 5 . 13 ) 
ov L 

In the accepted approximation in the rhs of this equation the difference between E and E' 
has to be neglected. 

According to our definition the energy E appearing at the Fourier transform of the 
transmission amplitude is the total energy of the system while E' is the energy characterizing 
the electron trajectory. The difference between them is the energy transferred to the phonon 
system, i.e., the average energy loss of the tunneling electron. 

The substitution of T = L/v and calculation of the integral with respect to q x in Eq . (|A8|) 
we obtain for the deformation potential, 

^dyn(E) = / |M q J 2 q± 2 d 2 q ± . (5.14) 

For the piezoelectric potential the terms in |H q „| 2 containing q x are small in a/L and the 
main contribution is 

* d yn(E) = / |2 q ,| 2 |M q± | 2 rfV • (5.15) 

For [100] tunneling direction Eq. ( |5.15| ) becomes 



oo2 2r2 „ „2„2 

*«nW = -^^TT- / |M q± | 2 ^ d\ ± , (5.16) 



and for [110] tunneling direction 



^vn(E) = 2 ~^r I l^xl 2 (i + 4g 2 g 2 2 ) ^ • (5.17) 
fju c j y I 
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It worthwhile to note that $<z yn is positive, i.e., it increases the transmission coefficient 
as well as <& sta t- One could expect that phonon emission makes the barrier for electron 
effectively higher that may lead to a reduction of the transmission probability. However, 
§dyn is calculated for E = E' and this effect can appear only in the next approximation. 

The average energy loss resulting from the deformation potential coupling is, 

E-E' = f^- / |M qx |V 2 dV . (5.18) 



87V 2 pV 2 o 

Whereas the average energy loss due to the piezoelectric coupling is given by: 



E-E' 



24/3VL f,^ |2 qlq 2 



/|M q J 2 MirfV, (5.19) 
J a i 



pv 2 e 2 J ' q~]_ 



for [100] tunneling direction and 

E - E ' = ^ / l^x! 2 (4 + ^ ■ (5.20) 

for [110] tunneling direction. 

The energy loss is proportional to the width of the barrier which means that it is accu- 
mulated along it. 

The comparison the static and dynamical phonon corrections to the tunneling probability 
gives &dyn/&stat ~ ( s / v )(L/a) where a is the width of the constriction. Typically the ratio 
L/a is not very large so the expansion that we used is justified. 



VI. DISCUSSION AND SUMMARY 

In this paper we have presented a detailed study of effect of coupling with acoustic 
phonons on electron tunneling across a rectangular barrier in a realistic situation. We stud- 
ied piezoelectric and deformation potential coupling at low temperature which means that 
Xt/cl ^> 1, where At is the thermal phonon wave length (see Appendix). We considered only 
a one-dimensional problem that can be realized in a quantum wire or a narrow constriction. 

In our calculation we assumed that the barrier is high enough to describe tunneling in 
the semiclassical approximation. Our detailed calculations reveal that the typical phonon 
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interacting with the tunneling electron is chosen through the length of the barrier L and the 
width of the constriction or quantum wire a. It is thus the geometry of the potential barrier 
which defines the phonon frequency u. Under such a condition the electron motion under the 
barrier is so fast that phonons don't follow it and can be considered nearly static during the 
time necessary for electron to traverse the barrier. The main effect of phonons in this case is 
a modulation of the barrier so that its height has to be considered as a quantum mechanical 
variable which probability distribution is controlled by a phonon wave function. In this 
case, roughly speaking, an electron chooses for tunneling the phonon configurations when 
the barrier is lower than its average value, that results in an increase of the transmission 
coefficient compared to that with zero electron-phonon coupling. Thus the interaction of an 
electron with the zero point phonon fluctuations increases the tunneling probability. 

The correction to the tunneling exponent in the static approximation is proportional to 
the length of the barrier. That is the coupling affects the dependence of the transmission 
coefficient on the height of the barrier only. 

The dynamical correction leads to two effects. First, it describes the reduction of the 
electron energy due to phonon emission. Second, which is probably more interesting, it 
gives the dependence of the transmission coefficient on the length of the barrier different 
from the regular one where the log of the transmission coefficient is linear in the length. We 
calculated only the first order correction to the exponential dependence. The effect can be 
stronger and is probably measurable for tunneling near the top of the barrier. One should 
note that even though the dynamical correction to the transition probability is smaller than 
the static correction, it is still an exponential correction and it easily can be made larger 
than unity, i.e., by changing the length of the barrier. In this case it can significantly affect 
the tunneling probability. 

The dependence of transmission coefficient on the height and the length of a barrier can 
be measured in devices where the barrier is introduced with the help of a gate. In such a 
device both the height and the length of the barrier are controlled by the gate voltage. Our 
results point out to a deviation of these dependencies from the standard ones obtained from 
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stationary Schrodinger equation. 

The application of our results to a long quantum wire can encounter a difficulty because 
we did not take into account electron-electron interaction. The effect of this interaction can 
be reduced for a short wire or a narrow constriction. 

We expect non-trivial results through further use of the piezoelectric and deformation 
potential coupling mechanisms in two or three dimensional systems. 

We believe that the physical considerations presented here which greatly simplified the 
calculations are convenient for extending the calculations to two and three dimensional 
physical situations. 
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APPENDIX: NON-STATIC CALCULATION OF THE TRANSITION 
AMPLITUDE AND TRANSITION PROBABILITY 



In the Appendix we carry out the integration with respect to initial and final phonon 
coordinates and expand the result in two small parameters defined in Sec.[TTT], s/v and 
Vint/iV — E). The calculation presented here are for the finite temperature case. 

The explicit form of Eq. ( |5.10| ) for T = T is 



£ ph {T-T) = g\T) 
x J duqi J duqj J duqi,exp 



q 



Z ^/ 



■PEr, 



tpm (Mqi)V'm 



qv ' 



where Z = J2m e ~ l3Em > ft is one over the Boltzmann constant times the temperature (we 
use this notation only in the Appendix and it should not be confused with the piezoelectric 
module in the body of the paper), 

Vq («q/, «qi, Uqi] T) = t q (| M qi f + | M q/ | 2 ) - (w qi W q/ + U^U^f 

+ e q (\ 1 2 + I Utf | 2 ) - (ucpU*^ + n qi w q/ 

+ uy* + U^I* + U~y 2 + Utfh - h , (A2) 



ot q = pu q /2smh(u> q T), e q = cosh(u q T), 
1 rL 



pU q V JO 
1 rL 



dxfq(x) sinh( 



0J q X s 



V 



pUJqV JO 

2 rL 



dxfq(x) sinh(^^ — 



dx dy fq{x)f*(y) + /*(z)/ q (y) sinh( ) sinh( 



p 2 ui q v 2 Jo Jo L q q 



Using the well known result for the density matrix |2l| : 



1pm (Uqt)V'm ( M q*) = eX P 



-ftl^PUql u w I tanh(^— , 



(A3a) 
(A3b) 
(A3c) 



(A4) 



where g((3) is given by: 



9((3) = U 



(A5) 



q y 2nh sinh hu q /3 

The result of the integration with respect to initial and final phonon coordinates gives 
i ph (T;T) = Aexp[2% h {T)/h] where 



a„ 



«0 = £? 



q e q 



e Q -h + h 



tq 2 + PSqjq/2 



+ 1 h 1 + 



(A6) 



7q = e g tanh(-^)/a 



Finite temperature effects in Eq . (|A6|) can be neglected if tanh(fru>j/3/2) pa 1. The physical 
meaning of this approximation is that the phonon thermal wave length is much larger than 
the quantum wires width i.e. Ay/a ^> 1. For the width of the constriction of 100 A or 
larger this condition is satisfied up to room temperature. From this estimate it is also follows 
that the matrix element for high energy phonons that could activate a tunneling electron to 
energy comparable with the height of the barrier is small and we neglect it. Thus the rest 
of the calculations will be performed for the zero temperature case. 

One should note that the phonon part of the transition probability, ^(TjT), is larger 
than 1, causing an enhancement of the electron tunneling probability. 

It is convenient to separate in 3> p h(T) a static part that is obtained by putting s/v = 
and a dynamical correction, 



%h(T) = $ stat (T) + $ dyn {T) . 



(A7) 



The static part 



stat 



E 



i 



dxfq(x) 



(A8) 



is equal to the expression appearing in Eq. (|4.7| ). The dynamical correction is necessary for 
the calculation of phonon emission and it is calculated in the leading order in s/v 1. For 
the deformation potential the dynamical correction is given by: 

sm(q x L) cos(q x L) 



4A 2 T 



q L 



9± 2 |MqJ 2 



(A9) 
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For the piezoelectric coupling one gets: 



$dyn{E) 



4T 



pv 2 V l 



vol 



sm(q x L) _ cos(q x L) 
q 3 x L q 2 x 



P I 2 IM I 2 

I 1 — 'qH | J "qjJ 



(A10) 
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